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I. ABSTRACT 

We demonstrate that full temporal characterisation of 
few-cycle electromagnetic pulses, including retrieval of 
the carrier envelope phase (CEP), can be directly ob¬ 
tained from Frequency Resolved Optical Gating (FROG) 
techniques in which the interference between non-linear 
frequency mixing processes is resolved. We derive a 
framework for this scheme, defined Real Domain-FROG 
(ReD-FROG), as applied to the cases of interference be¬ 
tween sum and difference frequency components and be¬ 
tween fundamental and sum/difference frequency com¬ 
ponents. A successful numerical demonstration of ReD- 
FROG as applied to the case of a self-referenced mea¬ 
surement is provided. A proof-of-principle experiment is 
performed in which the CEP of a single-cycle THz pulse 
is accurately obtained and demonstrates the possibility 
for THz detection beyond the bandwidth limitations of 
electro-optic sampling. 


II. INTRODUCTION 

Few-cycle electromagnetic pulses offer a means of gen¬ 
erating and controlling physical processes on sub¬ 
femtosecond timescales relevant in atomic and condensed 
matter physics. The characterisation of such short pulses 
remains a crucial aspect of the experimental process and 
has a rich history 1 3 . Of the many techniques available 
today, Frequency Resolved Optical Gating (FROG) has 
found widespread application due to its robustness and 
relatively simple experimental implementation. An ex¬ 
pansive array of related FROG techniques have been de¬ 
veloped in the last two decade^® each offer specific ad¬ 
vantages in varying experimental conditions, yet share 
a common mathematical and physical basis: the mea¬ 
surement of an intensity spectrogram from a non-linear 
combination of pulses. This is widely held to make de¬ 
termination of the carrier-envelope phase (CEP, (j ) CE ) 
beyond FROG techniqued^HUl 

For electromagnetic pulses with many cycles of the car¬ 
rier frequency the CEP approximates a simple time shift 
of the carrier oscillation and is often of limited physical 
interest; for few-cycle pulses however, CEP qualitatively 
alters the temporal profile of the electric field and subse¬ 
quently the interaction physics of the pulse. In current 


ultrafast metrology full temporal characterisation typi¬ 
cally requires a separate and often complex measurement 
of CEP in addition to information such as that provided 
by FROG. Stereo-ATI has been used to measure CEP by 
measuring the direction of photoelectrons emitted follow¬ 
ing ionisation in a time-of-flight spectrometer . 9 13 ^ CEP 
measurement has a lso be en demonstrated in so-called at- 
tosecond streaking j 16 * 17 * which utilises the spectral anal¬ 
ysis of photoelectrons excited by sub-femtosecond high 
harmonic pulses. In this latter case the principles of 
FROG have been applied to the emitted photoelectrons 
(as opposed to the input electromagnetic pulse as per 
usual) in the FROG-CRAB algorithm.® 

The CEP appears as the zero-order (constant) term in 
a Taylor expansion of phase. If <p CE is determined at one 
frequency it is therefore known uniquely. Recently, No¬ 
mura 10 and Shira 11 have utilised this, determining <p CE 
at low frequencies through THz electro-optic sampling 
(EOS) techniques, obtaining complete pulse information 
when combined with traditional FROG techniques. 

Here we show that self and cross-referenced FROG 
techniques can be extended to be directly capable of 
full characterisation of the electric field temporal pro¬ 
file of few-cycle electromagnetic pulses; similarly, we 
show and experimentally demonstrate the complete re¬ 
trieval of electric field can be obtained through cross- 
reference with higher frequency pulses, even when the 
higher-frequency probe is of similar duration and un¬ 
known carrier-envelope phase. There is no requirement 
for the to-be-determined pulse to contain low-frequency 
content accessible to ^-function like sampling, with re¬ 
trieval possible directly from frequency-mixing spectro¬ 
grams. 


III. RED-FROG: THEORY AND RETRIEVAL 
ALGORITHM 

FROG of electromagnetic pulse characterisation consists 
of a spectrally resolved auto- or cross-correlation of pulses 
mediated by a non-linear mixing process. Herein the de¬ 
tailed discussion assumes second-order mixing through 
X^ 2 \ although the approach and results can be readily 
modified to higher order mixing. For observable, real, 
electric fields in the time domain, Fi(t), F 2 (t) the com¬ 
plex spectra Fi(cc), F 2 (cc) are defined from the Fourier 
transform, introducing the mathematical construct of 
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negative frequencies. For the x^ interaction between 
pulses the spectrogram, the spectrum of x® generated 
field as a function of relative delay r between input 
pulses, can be given as: 


/(cj;t) = 



dQEi (cj 


2 


Q)E 2 (fi) exp (ifir) 


(1) 


in which we take the common approximation that the 
linear and non-linear material response functions can 
be collected outside the integral in the response func¬ 
tion R(uj )!^^ 9 \ Techniques such as surface harmonic 
generation 20 21 or phase-match angle-dithering 22 are ca¬ 
pable of very high bandwidth response, as would be nec¬ 
essary for few cycle pulses. As the specific form of R(uj) 
does not affect either analysis or conclusion it is omitted 
from the following discussion. 

The carrier envelope phase (j) CE is the phase contribu¬ 
tion that is widely accepted to be unmeasurable through 
FROG phase retrieval; to highlight conditions for obser¬ 
vation and retrieval, (p CE is explicitly separated out from 
the total spectral phase. All electric fields are explicitly 
constrained to be real in the time domain, which gives 


A'total (^-0 


E{uS) exp (icj) CE ) ; uj > 0 
E*(\uj\) exp(— i(j) CE ) ; uj < 0 


The imposition of the physical constraint for all fields to 
be purely real in the time domain and, in particular, the 
Heaviside functional form of (j) CE (uj ), is what gives rise 
to the observability of the carrier-envelope phase. Intro¬ 
ducing the functional form of equation [ 2 ] for the electric 
fields, the spectrogram can be expressed as 


/(w;r) = \SFG(lu;t)\ 2 + \DFG+(uj-,t)\ 2 + \DFG-(oj;t)\ 2 
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The contributions arising from sum-frequency ( SFG ) 
and difference frequency generation are represented by 
SFG and (DFG+, DFG- ) respectively. When both sum 
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FIG. 1. Schematic demonstrating bandwidth conditions for 
carrier envelope phase observation, for the two examples dis¬ 
cussed in the text. A) Self referenced mixing of a few cy¬ 
cle pulse, (top) spectrum for input pulse; (bottom) the sum 
and difference spectra. DFG+ = DFG- for this configura¬ 
tion. (B) Cross-correlation of spectrally distinct pulses, (top) 
input spectra; (bottom) frequency mixed output, for which 
DFG- =0. 


and difference frequency contributions spectrally overlap 
non-zero the absolute carrier phase becomes observable. 

In this work we address two specific cases; that of a 
single broad-band pulse cross-referenced with itself, and, 
as part of a proof-of-concept demonstration, the cross¬ 
correlation of two spectrally distinguishable pulses. The 
first case applies to the measurement of few-cycle optical 
pulses; the second case describes the unambiguous deter¬ 
mination of the electric field of a pulse through sampling 
with another (unknown) optical field, which may be of 
similar or longer duration than the lower frequency pulse. 
Representative field amplitude spectra of these two cases 
and the relative extent of the sum and difference fre¬ 
quency terms are shown schematically in Fig. [lj The 
principles discussed with these example can be readily ex¬ 
tended to address other specific spectral considerations, 
such as spectrally overlapping pulses. 

To demonstrate the complete electric field retrieval for 
the self-referenced measurement of a single pulse, test 
pulses with complex spectra differing only in carrier- 
envelope phase have been constructed, as shown in 


















3 



Time delay / fs 

-10 -5 o 5 10 



Time delay / fs 


FIG. 2. Numerical demonstration of the self-referencing retrieval of a transform-limited single-cycle optical pulse (110 THz 
bandwidth Gaussian with peak centred at 375 THz). A) Simulated experimental spectrogram with cj) CE = 0. B) Retrieved 
spectrogram, in which the retrieval has been performed over a truncated spectral region (frequencies above 375 THz); C) 
comparison of input and retrieved electric field profile. D-F) as above, with (j) CE = tt /2 rad. All spectrograms are shown with 
logarithmic intensity scale. 


Fig. |2p, F. Simulated ’experimental’ spectrograms were 
constructed directly as the product of the time domain 
fields. Modified FROG retrieval algorithms, accounting 
for the constraint of a strictly real time domain field, 
are then applied to the ’experimental’ spectrograms. Re¬ 
trieval has been examined in presence of both measure¬ 
ment noise and spectral truncation. 

Conventional FROG algorithms follow a two stage iter¬ 
ative procedure, the first being the comparison of the field 
amplitudes of the product spectrogram and the data. 
The second ’numerical-constraint’ draws on the physics of 
the relevant frequency mixing process to extract the indi¬ 
vidual time domain fields from the field-product spectro¬ 
gram. This constraint is conventionally undertaken from 
a functional minimisation between intensity spectrogram 
and the expected spectrogram inferred from iterated test 
fields and the physical process giving rise to frequency 
mixing. 

We have modified the PCGPA FROG retrieval algo¬ 
rithm 23 to include the physical constraint of a real time- 
domain field, which together with the frequency-mixing 
spectral overlap, we label as ReD-FROG (Real Domain- 
FROG) retrieval for distinction from conventional algo¬ 
rithms. The use of the PCGPA in providing the necessary 
deconvolution step decreases both the processing time 


and complexity of ReD-FROG, but is limited to measure¬ 
ments obeying equation [lj For more complex pulses, or 
instances for which dispersion of the mixing medium can¬ 
not be neglected, a more complex version of the decon¬ 
volution algorithm may be required.^ The spectral con¬ 
straint of real time-domain fields can be achieved through 
extending the measured spectrogram by mirror reflection 
about zero-frequency and applying the so-called ’numer¬ 
ical constraint’ minimisation procedure on the extended 
data set. Through effective inclusion of negative frequen¬ 
cies, and the pinning of the spectrogram to have mirror 
symmetry, the input field retrieved is constrained to be 
real in the time domain. In this case the mirror symme¬ 
try constrains the retrieved fields to have the Hermitian 
property necessary for the time domain fields to be purely 
real. 

The retrieval algorithm was applied on a spectral trun¬ 
cation of the measurement data; such a spectral limited 
data set is to be expected in any practicable implemen¬ 
tation. To fully investigate the effect of truncation, re¬ 
trieval was repeated for a range of cut-off frequencies be¬ 
tween 325 and 375 THz. 

As shown in Fig. |2p, F, the correct and full (including 
CEP) electric field temporal profile is recovered from the 
modified retrieval algorithm. All ReD-FROG retrievals 
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were based on a 512x512 grid. For Figs. HP and [2^, 
FROG errors of 5 x 10 -3 and 6 x 10 -3 were respectively 
obtained. In the given case the to-be-determined pulse 
was centred at 375 THz, and only the spectrogram data 
for frequencies > 375 THz was included in the retrieval; 
similar agreement in the retrieved profile was obtained 
across the full range of cut-off frequencies investigated. 
For an unambiguous electric field determination it is suf¬ 
ficient to include the sum-frequency region of the spec¬ 
trogram, in the presence of some overlapping difference- 
frequency mixing spectral content. Stable convergence 
was recorded with and without use of a spectral con¬ 
straint corresponding to a measurement of spectral in¬ 
tensity® applied over the truncated frequency range, al¬ 
though the rate of convergence was improved with its 
implementation. 

Up to this point in our analysis in we have taken 
the conventional form of the FROG spectrogram (Equa¬ 
tion [l]) , in which the only the non-linear generated fields 
are included. In this arrangement the bandwidth (A = 
c d u — uj l ) necessary for carrier envelope observability is 
A > 2uj l . By inclusion of the input field in the mea¬ 
surement process, additional (j) CE dependent interference 
terms arise which allow the bandwidth requirements to 
be relaxed. As a simplification, we assume that only one 
of the input fields spectrally overlaps with the mixing 
fields; under these conditions the spectrogram can be ex¬ 
pressed as: 


I(ui-,t)= 



2 


d QEi(oj — U)£ 1 2(^) exp(iUr) 


(5.) 


Explicit introduction of the CEP and expansion of the 
integral (as performed for equation |3| yields 


through interference between either the fundamental field 
Ei and the sum-frequency field, or between the funda¬ 
mental and Ei and difference frequency fields. The rela¬ 
tive bandwidth requirements are shown schematically in 
Fig. the carrier envelope phase is observable for a 
bandwidth satisfying A >uj l . This criterion is the same 
as that utilised in f-2f beating schemes for the CEP stabil¬ 
isation of repetitive pulse trains IH12M27 suc ] 1 stabilisation 
schemes however cannot provide a practical measurement 
of CEP, requiring extremely accurate characterisation of 
the dispersion associated with the experimental appara¬ 
tus. 3 14 Equation [6] demonstrates resolving this process 
as part of a FROG measurement provides the means to 
measure absolute phase without the need of a frequency 
comb. 

Figure [3^ presents an example of the spectrogram aris¬ 
ing from a self-referencing measurement of a transform- 
limited Gaussian pulse (75 THz bandwidth) in which 
overlap with the input field occurs. The overlap between 
sum and difference frequency generated fields is absent 
in this case and the interference features arise from the 
first three interference terms in equation [6] 
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+ 2 ft {sFG(u;t)DFG*_(lo-,t) e i2 ^ E } 



+25R | DFG+{oj] t)DFG*_ (w; t) e i2 ^ E ~ i2 4>2 E j 

The additional interference terms in the above expres¬ 
sion allow for carrier envelope phase observation even 
when there is no spectral overlap between sum and dif¬ 
ference frequency mixing. The first two interference 
terms in equation [6] yield measurement of the phase (j) EE 


FIG. 3(6)Schematic of bandwidth requirements for carrier 
phase observation through interference between fundamental 
and sum (or difference) frequency fields. (B,C) Spectrograms 
for a self-referencing measurement of a Gaussian transform- 
limited femtosecond pulse (375 THz carrier frequency, band¬ 
width 75 THz), <j) CE s= 0 and (j) CE = 7r/2 rad for (B) and (C) 
respectively. 
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IV. EXPERIMENT: CHARACTERISATION OF 
SINGLE-CYCLE THZ RADIATION 

We now consider the ReD-FROG process applied to the 
case of spectrally distinct pulses, as in Figure Eh For 
this configuration DFG-(uj\t) = 0, and therefore <$( E 
and <pi E — c p EE are not observable. Conversely the car¬ 
rier envelope phase (\> EE observable is independent of the 
(unknown) CEP phase of the upper band pulse. To ob¬ 
serve (j) EE we require Ai = 

In this case a benchmarked proof-of-concept experi¬ 
ment is feasible; a quasi-single cycle THz pulse can serve 
as the lower frequency pulse to be determined, while 
an optical pulse of similar bandwidth acts as the higher 
frequency probe. For the benchmarking, a significantly 
broader bandwidth ^-function like optical pulse can be 
employed for electro-optic sampling (EOS)pSEB obtain¬ 
ing an independent measure of the true electric field pro¬ 
file. 

Such an experimental demonstration has been per¬ 
formed, with the experimental arrangement shown in 
Fig.g An amplified femtosecond Ti:Sapphire laser sys¬ 
tem (Coherent Micra/Legend, amplifier output: 45 fs, 
lmJ at 800 nm) was used to drive a large-area semi- 
insulating GaAs photo-conductive antenna, producing 
quasi-single cycle sub-ps THz pulses. The pulses were 
detected within the near-field region of the antenna, al¬ 
lowing the electric field oscillation to have a unipolar ap¬ 
pearance. The optical probe pulse was selected prior to 
the PC A using a beam-splitter and directed through a 
zero-dispersion 4-f filter, allowing the pulse bandwidth to 
be selected as required. The benchmarking measurement 
of the temporal profile of the THz pulse was performed 
using a standard balanced detection EOS arrangement 
using a 45 fs optical probeP^ For the spectrogram re¬ 
trieval experiments the optical bandwidth was restricted 
to 1 THz, less than the full bandwidth of the THz pulse, 
with corresponding increase in duration to 500 fs. Optical 
probe temporal envelopes were confirmed with a custom- 
built SHG-FROG setup with 100 /im thick BBO. 

Both optical and THz beams were focused into a 
(110) orientated ZnTe crystal. The frequency-mixing sig¬ 
nal was isolated from residual input-probe light using 
a Gian-Thompson calcite polariser and spectrally anal¬ 
ysed as a function of the relative delay between the THz 
and optical probe beams using an imaging spectrometer 
(Jobin Yvon, iHR550) and intensified CCD camera (PCO 
DiCamPRO). A ReD-FROG algorithm as described ear¬ 
lier was applied to retrieve the field, using only the mea¬ 
sured spectrogram as input. 

The measured and retrieved spectrograms are shown 
in Figures [4)3, C, together with the ReD-FROG retrieved 
field time profile. The FROG error in the retrieved spec¬ 
trogram was 0.01. Despite the complex form of the THz 
pulse, both main peak and oscillations over a period 
of several picoseconds (due to THz absorption by am¬ 
bient water vapour) are recovered almost perfectly and 



Time delay / ps Time delay / ps 

FIG. 4. A) Schematic of experimental arrangement for car¬ 
rier phase retrieval and benchmarking. B) Measured spectro¬ 
gram with a 1 THz bandwidth, 500 fs duration, optical probe; 
C) ReD-FROG retrieved spectrogram D) comparison of true 
electric field as determined by electro-optic sampling with a 
45 fs probe, and that from the ReD-FROG retrieval. The 
electro-optic sampling measurement for the spectrally nar¬ 
rowed 500 fs probe is also shown for comparison. Spectro¬ 
grams are shown with logarithmic intensity scale. 

the CEP of the single-cycle pulse is correctly identified. 
For comparison the field profile inferred by direct electro¬ 
optic sampling with the 1 THz bandwidth probe, without 
spectrally resolved phase retrieval, is also shown. 

Whilst this experiment confirms the validity of CEP 
retrieval using ReD-FROG, the approach also finds im¬ 
mediate application in surpassing probe-duration limi¬ 
tations present in high time-resolution EOS. Continuing 
the notation of section II, these limitations can be ex¬ 
pressed in the frequency domain as A uj h ; for the same 
optical probe duration, we estimate that ReD-FROG 
(A > uj l ) is capable of measuring THz pulse durations 
x3-4 shorter than that accessible through EOS. The only 
existing or proposed THz schemes with similar capabil¬ 
ities, such as spectral upconversion 31 or BMX-FROGP 
rely on multi-stage detection in which the output of the 
electro-optic effect is measured as part of a second, dis¬ 
tinct non-linear process (such as an auto-correlation or 
FROG). The use of consecutive non-linear stages for 
pulse analysis limits the available signal to noise ratio and 
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thus the range of THz sources that can be measured; the 
application of the ReD-FROG method, being a single- 
stage direct analysis of the EO effect, overcomes such 
limitations and presents as a more robust and versatile 
measurement technique. 


V. CONCLUSIONS 

We have demonstrated that a spectrally-resolved mea¬ 
surement of the superposition of non-linear frequency 
mixing processes is capable of unambiguous electric field 
temporal characterisation, including the measurement of 


carrier envelope phase. Numerical simulations have been 
used to demonstrate the capability for few-cycle optical 
pulses. A proof-of-concept benchmarked experiment us¬ 
ing single-cycle THz radiation was performed and con¬ 
firmed the validity of the formalism. This method opens 
up new avenues of possibility in the measurement of few- 
cycle electric fields and demonstrates that, contrary to 
previous expectations, FROG is directly capable of mea¬ 
suring absolute phase. Potential applications for ReD- 
FROG in high time-resolution electro-optic sampling are 
also envisaged as a means to exceed the bandwidth lim¬ 
itations of standard electro-optic sampling techniques. 
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